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ABSTRACT 



Aims. The stability of dissipative Taylor-Couette flows with an axial stable density stratification and a prescribed azimuthal magnetic 
field is considered. 

Methods. Global nonaxisymmetric solutions of the linearized MHD equations with toroidal magnetic field, axial density stratification 
and differential rotation are found for both insulating and conducting cylinder walls. 

Results. Flat rotation laws such as the quasi-Kepler law are unstable against the nonaxisymmetric stratorotational instability (SRI). 
The influence of a current-free toroidal magnetic field depends on the magnetic Prandtl number Pm: SRI is supported by Pm > 1 
and it is suppressed by Pm J 1. For too flat rotation laws a smooth transition exists to the instability which the toroidal magnetic field 
produces in combination with the differential rotation. This nonaxisymmetric azimuthal magnetorotational instability (AMRI) has 
been computed under the presence of an axial density gradient. 

If the magnetic field between the cylinders is not current-free then also the Tayler instability occurs and the transition from the 
hydrodynamic SRI to the magnetic Tayler instability proves to be rather complex. Most spectacular is the 'ballooning' of the stability 
domain by the density stratification: already a rather small rotation stabilizes magnetic fields against the Tayler instability. 
An azimuthal component of the resulting electromotive force only exists for density-stratified flows. The related alpha-effect for 
magnetic SRI of Kepler rotation appears to be positive for negative dp/dz < 0. 
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1. Motivation 

This work is motivated by theoretical and experimen- 
tal progresses in studies of the stratorotational instability 
(SRI) and the magnetorotational instability (MRI) in MHD 
Taylor-Couette experiments. It has been shown theoretically 
(Molemaker, Mc Williams & Yavneh 2001 ; Yavneh, Mc Williams 
& Molemaker 2001; Dubrulle et al. 2005; Shalybkov & Rudiger 
2005; Umurhan 2006) and in the laboratory (Le Bars & Le Gal 
2007) that a combination of a centrifugal-stable nonuniform ro- 
tation law and a stable axial density stratification leads to the 
so-called stratorotational instability (SRI) in the Taylor-Couette 
flow. This instability exists only for nonaxisymmetric distur- 
bances. On the other hand, there are also nonaxisymmetric insta- 
bilities for a combination of Rayleigh-stable rotation laws and 
azimuthal magnetic fields (Pitts & Tayler 1985). The question 
is whether the combination of density stratification, differential 
rotation and toroidal fields acts stabilizing or destabilizing or 
whether even new instabilities arise. 

Such a combination of axial density stratification, stable ro- 
tation law and strong toroidal magnetic field is the typical con- 
stellation in accretion disks (Ogilvie & Pringle 1996; Curry & 
Pudritz 1996; Papaloizou & Terquem 1997). There the rota- 
tion is Keplerian with Q oc R^l 2 and the toroidal field is gen- 
erated from weak large-scale poloidal fields by induction due 
to the differential rotation. The standard case is that the result- 
ing toroidal field strongly exceeds the amplitude of the original 
poloidal field if the magnetic Reynolds number of the differen- 
tial rotation Rm = QR 1 jr\ is much larger than unity. If this is true 
then the standard MRI, i.e. the influence of the poloidal field on 
the stability of the differential rotation, would be of minor im- 



portance The question is whether the toroidal magnetic field 
can reach such high amplitudes or whether it becomes unstable 
already for much smaller values. It is known that toroidal fields 
with strong electric currents become unstable (Tayler instability, 
'TP ) but it is also known that in combination with differential ro- 
tation also toroidal fields become unstable which are current-free 
in the considered domain (azimuthal MRI, 'AMRI', see Rudiger 
et al. 2007a). In the latter case, with oc /J -1 , the questions 
are whether the density stratification destabilizes AMRI and/or 
whether the toroidal field stabilizes the SRI too strongly so that 
its real existence becomes basically suppressed. 

We shall consider the interaction of the differential rota- 
tion with both an axial density stratification and a toroidal mag- 
netic field in the simplified Taylor-Couette geometry. The den- 
sity stratification is always supposed to be stable but the mag- 
netic field and the rotation law between the cylinders can be both 
stable or unstable. If, in particular, the toroidal field is Tayler un- 
stable then the interaction of differential rotation, density stratifi- 
cation and magnetic field becomes highly complex. We find sta- 
bilization and destabilization in strong dependence on the mag- 
netic Prandtl number. Again, in experiments, for small magnetic 
Prandtl number, the flows are predicted to be stabilized. For 
galaxies and protoneutron stars (PNS) with their high magnetic 
Prandtl numbers we find the opposite: the magnetic influence 
supports the SRI leading to even smaller Reynolds numbers than 
in hydrodynamics. 



1 the same is true for calculations about the stability of strong 
poloidal fields without or with rotation 
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Fig. 1. The geometry of the problem - two concentric cylinders 
with radii R m and R out rotating with Q m and Q out . B<p is the az- 
imuthal magnetic field which, generally, is produced by both 
an axial current inside the inner cylinder and an axial current 
through the fluid. 

2. The Taylor-Couette geometry 

A Taylor-Couette container is considered which confines a 
toroidal magnetic field with given amplitudes at the cylinder 
walls which rotate with different rotation rates Q (see Fig.[TJ. 

The fluid between the cylinders is assumed to be incompress- 
ible and dissipative with the kinematic viscosity v and the mag- 
netic diffusivity 77. Derived from the conservation of angular mo- 
mentum the rotation law Q(R) in the fluid is 



Q = a H — r 

R 2 

with 



a = — TTT^ir 



1 -fj 2 



and 
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where cylindric coordinates (R,(p,z) are used and i2; n and Q out 
are the imposed rotation rates at the inner and outer cylinders of 
radii R[ n and R out . 

Similarly, the magnetic profiles are restricted by the az- 
imuthal component of the induction equation to 



Bi, - AR + 



B 

R' 



(4) 



where the first term corresponds to a uniform axial electric cur- 
rent at radius R and the second term is current free for R > 0. In 
analogy with /jq it is useful to define the quantity 



Mb 



B i 

~b~, 



AR out + BR ^ t 



AR;„ + BR 



1 ' 



(5) 



measuring the variation in B^ across the gap between the cylin- 
ders. The coefficients A and B are 
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The MHD equations for incompressible stratified fluids are 
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div 11 



+ (uV)p = 0, 



0. 



div B = 0, 



(7) 



where u is the velocity, B is the magnetic field, P is the pressure, 
g is the gravitational acceleration (supposed as vertical and con- 
stant), v is the kinematic viscosity, 77 is the magnetic diffusivity 
and yUo the magnetic constant. 

In the presence of a vertical density gradient (p = p(z)) it has 
been shown for B = 0, that the system (O allows the angular 
velocity profile ([TJ only in the limit of slow rotation and small 
stratification 
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It easy to show that the same is true for B + 0. So, we are inter- 
ested in the stability of the basic state 



U = (Q,RQ(R),0), 
P = Po(R) + P 1 (R,z), 



B = (0,B,(tf),0), 
p = po +Pi(z), 



(9) 



where Q is given by (Q}, B^ by (O, po is the uniform reference 
density, P is the total pressure including the magnetic part with 
\P l /P Q \ « 1 and| Pl /p | « 1. 

The linear stability problem is considered for the perturbed 
state of U +u, B+b, po+pi +p', Pq+Pi+P' . Using the conditions 
(O the linearized system (0 takes the Boussinesq form with the 
coefficients depending only on the radial coordinate, and a nor- 
mal mode expansion of the solution F = F{R)exp{i(m<p+kz+u>t)) 
can be used, where F represents any of the disturbed quantities. 

Finally after a normalization we find 
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where the same symbols are used for the normalized quanti- 
ties except P' which denotes P'/po and redefining p' as PmRe 
gp'/po- The dimensionless numbers of the problem are the 
magnetic Prandtl number Pm, the Hartmann number Ha, the 
Reynolds number Re 



Pm = -, 



Ha = 



B;„R() 



and the buoyancy frequency N after 



Re = 



8 dpi 
Po dz 



(11) 



(12) 



We used Po = (Pin (Pout - Pin)) 1 ' 2 as the unit of length, rj/Ro as 
the unit of the perturbation velocity, B[ n as the magnetic field 
unit, Q m as the unit of u, N, and Q, RqQ 2 as the unit of g, po as 
the density unit and vrj/R^ as the unit of the redefined P' , It is 
convenient to describe the influences of the density stratification 
and the basic rotation by the Froude number 



„ ^in 

Fr = — . 

N 



(13) 



A detailed description of the numerical methods has been 
given in earlier papers (see e.g. Riidiger et al. 2007b) and will not 
be reproduced here. Always no-slip boundary conditions for the 
velocity on the walls are used. The tangential electrical currents 
and the radial component of the magnetic field vanish on the 
conducting walls. For insulating walls the magnetic field must 
match the external magnetic field. 



3. Hydrodynamics 

According to the Rayleigh condition the Taylor-Couette flow 
with the rotation law (Q~|i is stable for 



(14) 



This condition has been extended to stratified fluids by Ooyama 
(1966). The general condition for nonaxisymmetric solutions is 
not known. It is known, however, that the axisymmetric mode is 
the most unstable mode for dissipative Taylor-Couette flows for 
pn > and the nonaxisymmetric modes can be more unstable 
than the axisymmetric ones only for counter-rotating cylinders 
(Langford et al. 1988). 

Surprisingly, a stable vertical density stratification destabi- 
lizes the Taylor-Couette flow with decreasing angular velocity 
(fin < 1) even beyond the Rayleigh line. This stratorotational 
instability (SRI) is nonaxisymmetric, the most unstable mode is 
the mode with m = I. Shalybkov & Riidiger (2005) found that 
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Fig. 2. The marginal stability line for m — 1 disturbances in a 
container with fj = 0.5 and /jq = 0.3. The instability is stabilized 
for both too weak and too strong stratifications. The minimum 
Reynolds number belongs to a Froude number of 1 .4 depending 
on the rotation law. 



the instability condition is /in < 77 instead of /in < 1 which 
has experimentally been confirmed by Le Bars & Le Gal (2007). 
This finding based on a restricted number of calculations. More 
detailed calculations are here reported with new results extend- 
ing the previous conclusions. 

In Fig. [2] the marginal stability line (the line which separates 
stable and unstable regions) is calculated as function of N 2 /Q 2 n 
for m = 1 disturbances and for a standard container with f) = 0.5 
and with /in = 0.3. The latter value lies beyond the Rayleigh line 
where nonstratified flows are stable. 

We are looking for the minimum critical Reynolds number 
which exists for a buoyancy frequency NIQ- m 0.71, or - what 
is the same - Fr 1 .4. Figure [2] reveals the SRI as a delicate 
balance of buoyancy and rotation law. Both the frequencies must 
approximately be equal otherwise one of the two stabilities dom- 
inates the constellation. 

The question is how this balance depends on the rotation 
law and on the geometry of the container. In Fig. [3] three con- 
tainers are used with different gap sizes. The computations also 
concern various density stratifications. For any gap size three 
vertical lines are given. Left is the line for the Rayleigh limit 
where Q oM = i3; n P 2 n /P 2 ut . The right line mimics the rotation 
law Q oc l/R which is typical for galaxies (i2 out = £?i n Pin/Pout)- 
The central line represents the Kepler rotation law Q oc P -1 - 5 , i.e. 
//n = ff , which we shall call the pseudo-Kepler limit. The main 
result of the Fig. [3] is that the pseudo-Kepler rotation is always 
unstable. This, however, is not true for the more flat pseudo- 
galactic rotation which for large gaps becomes stable. From the 
bottom plot in Fig.[3]we have to realize that the pseudo-galactic 
rotation law with Q oc 1 /P is already too flat for the possible 
existence of nonmagnetic SRI in galaxies. 

The present calculations do not confirm the p — f) line as 
the limit for the SRI. This limit can be larger for small-gap con- 
tainers and can be smaller for wide-gap containers. The stability 
limit seems to approach the Rayleigh line in the limit of very 
wide gaps and goes to unity for very narrow gaps. We shall ask 
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in the following whether such results are influenced by the exis- 
tence of a toroidal magnetic field. 



4. SRI with azimuthal magnetic field 

Without density stratification the ideal Taylor-Couette flow with 
imposed azimuthal magnetic field is stable against axisymmetric 
disturbances if 
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Fig. 3. The marginal stability lines for small-gap (top, fj = 0.78), 
modest-gap (fj = 0.5, middle) and wide-gap (fj = 0.3, bottom). 
The vertical dotted lines denote the Rayleigh limit, the pseudo- 
Kepler rotation and the pseudo-galactic rotation (fin = fj)- The 
curves represent different density stratifications with Fr = 2.2 
(dot-dashed), Fr = 1 (dashed), Fr = 0.7 (dotted), Fr = 0.5 (solid, 
maximum stratification). 
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(Michael 1954). This condition is the combination of the 
Rayleigh stability condition for differential rotation without 
magnetic field and the magnetic field stability condition with- 
out rotation. Without a global rotation the magnetic field is sta- 
ble (the second term in (15[ is positive) for the magnetic field 
profile © if 

< hb < - = £0. (16) 

Therefore, all magnetic profiles with positive fig or with fig < 2 
(for fj - 0.5) are stable against axisymmetric perturbations. As 
the dissipative effects are stabilizing the flow a magnetic field 
with jig beyond the interval dT6b becomes unstable against ax- 
isymmetric disturbances in case that the magnetic field ampli- 
tude (or its Hartmann number) exceeds a critical value. 

The stability condition for the azimuthal magnetic field 
against nonaxisymmetric disturbances after Tayler (1961) is 



dR 



(RBl) < 0. 



For the magnetic field profile dU it takes the form 



<Hb < 



4t)(1 - /) 2 ) 
3 - 2fj 2 - ff 



(17) 



(18) 



for the m — 1 (kink) mode which is the most unstable mode. 
Hence, all azimuthal magnetic fields with positive fig smaller 
than 0.62 (for fj = 0.5) are stable against nonaxisymmetric dis- 
turbances. Current-free fields with B^ cc l/R, i.e. jiB = 0.5 for 
fj = 0.5 are thus always stable against m — and m = 1. 

As always p.\ < fiQ the stability interval ( fTSI l is smaller than 
the stability interval ( fT~6b . In this sense the nonaxisymmetric 
('kink') disturbances are more unstable than the axisymmetric 
('sausage') disturbances. However, the situation is much more 
complex under the presence of (differential) rotation and which 
instability really dominates depends on the parameters of the 
problem (Riidiger et al. 2007a,b). 

In the present paper the influence of stable vertical density 
stratifications on the Taylor-Couette flow stability with imposed 
azimuthal magnetic fields is considered for three different cases: 

1) both the rotation law and the magnetic field are individually 
stable and both together are unstable (i.e. hb = 0.5, no criti- 
cal Hartmann number) 

2) the magnetic field is stable for m = but unstable for m = 1 
(i.e. /is = 1, one critical Hartmann number) 

3) the magnetic field is so steep that it is unstable for m-Q and 
for m — 1 (i.e. jig = 3, two critical Hartmann numbers). 



4.1. AMRI with density stratification 

The azimuthal magnetic field with cc l/R (i.e. jis = 0.5 
for fj = 0.5) is considered which is current-free between the 
cylinders as the simplest example of an azimuthal magnetic field 
which is basically stable without rotation for both axisymmetric 
and asymmetric disturbances. 
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The typical stability diagram is presented by the Figs.|4]and|5] 
for containers with conducting/isolating cylinders, medium gap 
width (fj = 0.5) and rather fiat rotation law (pseudo-Kepler ro- 
tation with fiQ - 0.35 and a slightly steeper rotation law with 
fi = 0.3) which are stable without density stratification and with- 
out magnetic field. We find that the results depend only slightly 
on the conducting properties of the cylinder material. 
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Fig. 4. The marginal stability lines for the flow with conducting 
cylinders, f) = 0.5, current-free (p B = 0.5) azimuthal magnetic 
field, pseudo-Kepler rotation (pn = 0.35), Fr = 0.5 for distur- 
bances with m = 1 (solid) and m = 2 (dotted). The curves are 
labeled by their magnetic Prandtl numbers. For comparison the 
dashed line demonstrates AMRI (m =1) for homogeneous fluids 
(N = 0). 



The solid line in Figs. [4] and [5] is the marginal stability line 
of a homogeneous fluid with Pm = 1 . There are no solutions at 
both the vertical axis and the horizontal axis. The MHD flow is 
only unstable as a combination of differential rotation and the 
azimuthal magnetic field. For any supercritical Hartmann num- 
ber there are two critical Reynolds numbers between which the 
fluid is unstable. Similarly, for any supercritical Reynolds num- 
ber there are two critical Hartmann numbers between which 
the fluid is unstable. We have called this phenomenon as the 
Azimuthal MagnetoRotational Instability (AMRI) which scales 
with the magnetic Reynolds number in the same way as it is 
known from the standard MRI in Taylor-Couette experiments. 
However, magnetic Reynolds numbers of order 100 are too high 
to be realized in the MHD laboratory. 

With a density stratification (Fr = 0.5) the situation changes. 
The following results can be interpreted as either AMRI for 
density-stratified fluids or as the influence of an azimuthal mag- 
netic field on SRI. There is now an instability (SRI) without any 
magnetic field for Re > 296 at m - 1 and for Re > 366 at m = 2 
(for conducting cylinders, pseudo-Kepler rotation law, Fig. |4]i 
and for Re > 264 at m = 1 (insulating cylinders, pa = 0.3, Fig. 
[5). For this latter case the calculations are restricted to the kink 
(m = 1 ) instability. 

For conducting boundaries (see Fig. |4j and for a special case 
(Ha = 100, m = 1) the radial eigenfunctions are given for 
Pm = 1 (Pig. [6J and Pm =10 (Pig. |7j which demonstrate that 
the SRI modes do not form boundary phenomena. The profiles 
are normalized since in the linear theory the amplitudes have no 
own physical meaning. For the linear growth time in units of the 
rotation period of the inner cylinder for the solution with Pm — 1 
one obtains 



^"growth 



243 



Re - 341 



(19) 



Hence, for Re = 400 the e-folding time is about four rota- 
tions. The vertical wave number at the critical Reynolds number 
Re = 341 is 8.63 what means that the vertical cell size in units 
of the gap width is about 0.36. Without magnetic field the wave 
number is 8.16. Common action of magnetic field and density 
stratification leads to rather flat cell configurations. 
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Fig. 5. The same as Fig.|4]but for insulating cylinders with pn 
0.3 (as in Fig. [2j but without m = 2 mode. 
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Fig. 6. The eigenfunctions for the radial components of flow 
(left) and field (right) for Pm = 1 . Solid: real part, dashed: imag- 
inary part; the (arbitrary) amplitude is normalized. Fr = 0.5, 
Ha= 100, m= 1, see Fig. H 



Note that the influence of the azimuthal magnetic field on 
the SRI strongly depends on the magnetic Prandtl number Pm. 
For Pm < 1 the SRI is stabilized by the magnetic field but it 
can be supported for Pm > 1 (see Figs. [4] [5}. For Pm = 10 
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Fig. 7. The same as in Fig. [6]but for Pm = 10. The eigenfunc- 
tions do not reflect a boundary phenomenon. 



H 1000 



the critical Reynolds number is reduced (by a factor of two) for 
Ha a; 200. In this case AMRI supports the SRI. This behavior 
is similar to that of the standard MRJ with axial magnetic field 
and hydrodynamically unstable rotation where also for Pm > 1 
the magnetic field supports the centrifugal instability (Riidiger 
& Shalybkov 2002). 

Note also that the magnetic suppression of the SRI for small 
Pm is rather weak. Even for Pm = 10~ 5 a (large) Hartmann num- 
ber of order of 100 only leads to a small critical Reynolds num- 
ber of about 500. The lines for Pm = 10~ 5 in Figs. [4] and [5] are 
identical with the lines for all smaller Pm, i.e. they also hold 
for gallium (10~ 6 ) and also mercury (10~ 7 ). Here all the curves 
scale with Re rather than Rm so that the corresponding mag- 
netic Reynolds numbers are quite small compared with those 
for AMRI. As also the corresponding Hartmann numbers with 
Ha 5 100 are not too high, laboratory experiments with fluid 
metals should be possible provided a sufficiently density strat- 
ification can be produced. 

It is also shown in Fig. [4] that for a given Pm the critical 
Reynolds numbers for m — 2 exceed those for m = I. The differ- 
ences, however, are not too big so that for slightly supercritical 
Re several modes should be excited. This is insofar surprising as 
the known smoothing action of differential rotation upon modes 
with high m here seems to be very weak. 

The rather small Reynolds numbers and Hartmann numbers 
lead to the impression that the instability should be observable 
into the laboratory. For hydrodynamically unstable situations we 
have only a suppression of the instability by the magnetic field 
for small Pm. For hydrodynamically stable situations (solid lines 
in Figs. |4] and |5]l it was shown for uniform fluids that the appro- 
priate numbers of the problem are the magnetic Reynolds num- 
ber Rm and the Lundquist number S 
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Fig. 8. AMRI for stratified fluids with Fr = 0.5 between insu- 
lating cylinders (77 = 0.5, hb = 0.5). The rotation law is so flat 
(yUn = 0.6) that nonmagnetic SRI does not exist (see Fig. [3] mid- 
dle). Lines are labeled by their Pm numbers. 



The solid line in Fig. [9] demonstrates that the AMRI (like 
the standard MRI) exists for all rotation laws with decreasing 
angular velocity as a function of radius (i.e. for jun < 1). 
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rather than Re and Ha. Figure [8] demonstrates that the same is 
true for density-stratified fluids with so flat rotation laws that hy- 
dromagnetic SRI does not operate (see Fig.0. Both the critical 
magnetic Reynolds number and Lundquist number do not vary 
remarkably when the magnetic Prandtl number varies over 3 or- 
ders of magnitude. 

Stratified flows without magnetic field become stable for /jq- 
values slightly greater than /in = 0.5. According to Fig. [8] the 
flow is stable for /in = 0.6 without a magnetic field but it be- 
comes unstable with magnetic fields. The critical Rm are not 
monotonous with Pm. The values of Rm decrease with Pm for 
small Pm and v. v. The result can be seen in Fig. [8] which differs 
in this respect to the AMRI without density stratification (see 
Riidiger et al. 2007a). 
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Fig. 9. The minimum Reynolds numbers for insulating cylin- 
ders, 77 = 0.5, Pm = 1 and for the kink (m = 1) mode. The verti- 
cal lines are the same as in Fig. [3] The numbers on the curves are 
Hartmann numbers which correspond to the minimum Reynolds 
numbers. The solid line is AMRI without a stratification, the dot- 
dashed line is AMRI with stratification (Fr = 0.5, see Fig. [8]) and 
the dashed line is SRI without magnetic field. 



According to the Figs. [4] and [5] the marginal stability lines 
have always a minimum for some Hartmann number including 
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Ha = 0. These minimum Reynolds numbers are plotted in Fig. [9] 
as functions of fin. The dashed line is SRI without magnetic field 
which disappears for too flat rotation laws (yUn > 0.52). There is 
a smooth transition, however, to the instability AMRI (current- 
free toroidal field plus differential rotation) which has been com- 
puted here for the first time with a density gradient (dot-dashed). 
The comparison with standard AMRI without dp/dz (solid line) 
yields the expected suppression by density stratification. Note, 
however, that for quasi-galactic rotation (fici = 0.5 in Fig.|9]l the 
SRI with magnetic field can easier be excited than without mag- 
netic field. As the figure shows the necessary magnetic field has 
a Hartmann number of about 137. 

As AMRI scales with Rm rather than Re. the 'smooth' tran- 
sition interval from the hydrodynamic solutions to the magneto- 
hydrodynamic scales with Pm -1 . Hence, the dot-dashed line in 
Fig.[9]becomes more and more steep for decreasing Pm. 

4.2. Density stratification and Tayler instability (m = 1) 

The influence of the vertical density stratification on the non- 
axisymmetric Tayler instability (i.e. fi\ < fig < P-o) is consid- 
ered for rotation laws which are steep enough to be Rayleigh 
unstable. Again the gap width is fj = 0.5, the outer cylinder is 
resting (fiQ = 0) and the toroidal field is rather homogeneous 
(jj-b = 1) so that without rotation only the mode m = 1 is unsta- 
ble. Calculations have been made for both insulating and con- 
ducting cylinders (Fig.fTOb. The difference between conducting 
and insulating boundary conditions proves to be only quantita- 
tive. 

For m — 1 the critical Hartmann numbers above which 
the flow is unstable for Re = are 150 for conducting cylin- 
ders and 109 for insulating cylinders. These critical numbers 
are not influenced by the density stratification. For m — no 
critical Hartmann numbers here exist; the magnetic field stabi- 
lizes this Rayleigh-mode for both homogeneous (dotted line) and 
density-stratified (dot-dot-dot-dashed lines) fluids. For homoge- 
neous fluids the (dotted) marginal stability line for axisymmetric 
disturbances does not depend on the magnetic Prandtl number 
Pm (Shalybkov 2006)! 

The critical Reynolds numbers for the Rayleigh instability 
for Ha = do not depend on the conducting properties of the 
cylinders. They are 68 (m = 0) and 75 (m = 1) without stratifi- 
cation. The density-stratification stabilizes the flow; the critical 
Reynolds numbers are 294 for m = and 226 for m = 1 (for 
Fr = 0.5). With density-stratification the kink mode (m = 1) 
proves to be the most unstable mode. 

Moreover, for given Pm the kink mode is the most unstable 
one for all Ha for Fr = 0.5. For small Pm this mode is stabi- 
lized by weak magnetic fields before a dramatic destabilization 
happens for larger magnetic fields. For increasing Pm the flow 
becomes unstable for decreasing Reynolds numbers at a fixed 
Hartmann number. The opposite is also true: the smaller the Pm 
the stronger is the magnetic field influence and the flow becomes 
unstable for smaller Hartmann numbers at a fixed Reynolds 
number. 

Note that the competition between magnetic and centrifu- 
gal instabilities can lead to a rather complex behavior as illus- 
trated by the marginal stability line with Pm = 10. Generally, 
however, the stability region is much larger for density-stratified 
fluids so that we find that a stable stratification increases the 
flow stability. Of particular interest is the phenomenon that for 
density-stratified fluids generally the critical Hartmann numbers 
are larger than without rotation and the same is true for small Pm 
with respect to the critical Reynolds number ('ballooning'). For 
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Fig. 10. The marginal stability lines for a flow with insulating 
(top) and conducting (bottom) cylinders. It is 77 = 0.5, fin = 
(resting outer cylinder) and fig = 1. The curves are labeled 
by their magnetic Prandtl Pm. The m = (dotted) mode (no 
dependence on Pm) and m = 1 (solid) mode for homogeneous 
fluid and m = (dot-dot-dot-dashed) and m = 1 (dot-dashed) 
modes for density-stratified fluids with Fr = 0.5. 



large magnetic Prandtl number, however, the magnetic field ba- 
sically destabilizes the Rayleigh instability. Even for weak fields 
the critical Reynolds number is smaller than the Reynolds num- 
ber for the nonmagnetic case. 

4.3. Density stratification and Tayler instability (m = Q,l) 

The most complex situation appears when both the modes with 
m — and m = 1 are Tayler-unstable. Only conducting cylin- 
ders with medium-sized gap (fj = 0.5), outer cylinder at rest 
and a toroidal magnetic field with strong currents (jig = 3) are 
analyzed. The differential rotation is Rayleigh unstable without 
magnetic field and density stratification and the magnetic field is 
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Fig. 11. The marginal stability lines for conducting cylinders 
with f) = 0.5, /J.n -Q,fi B - 3, Pm = 1(T 5 (top) and Pm = 1 (bot- 
tom). The lines connect the Rayleigh instability without mag- 
netic field with the Tayler instability without rotation. The dotted 
lines are for m = 0; the solid lines for m — 1 . The inner curves 
are for homogeneous fluids (N = 0) and the outer curves for 
stratified fluids with Fr = 1. The density stratification stabilizes 
in both directions. 



Tayler-unstable without rotation (for m — and m = 1). The crit- 
ical Hartmann number for Re = is nearly equal for all cases: 
m = gives Hao = 21 and m = 1 gives Hai = 20.6. For other pa- 
rameters it can be more different. Either the m = or the m = 1 
mode can be the most unstable mode. Often also the m = mode 
is the most unstable for small Ha numbers and fast differential 
rotation (see Riidiger et al. 2007b). 

The critical Reynolds numbers above which the rotation be- 
comes unstable without a magnetic field (Ha = 0) are called Reo 
and Rei . Reo is smaller than Rei for flows with < fin < ff. The 
density stratification increases Reo and Rei (Fig.fTTI). Note that 
with density-stratification the m = 1 mode becomes the most 
unstable mode. With magnetic field it is also the most unstable 
mode almost everywhere except for Pm = 10~ 5 and Ha > Hao. 

Note that 'ballooning' of the stability region is produced by 
the density stratification. The density stratification in combina- 
tion with the basic rotation stabilizes the flow for larger mag- 
netic fields. Under the influence of differential rotation and den- 
sity stratification stronger magnetic field amplitudes prove to be 
stable than they were without rotation. The effect already exists 
for rather slow rotation rates (see Fig. [TT] Pm =1). This phe- 
nomenon is a consequence of the explicit inclusion of the sta- 



ble density stratification into the MHD equations. One finds that 
both the Rayleigh instability and the Tayler instability are stabi- 
lized by the density stratification. Without rotation, however, the 
density influence of the Tayler instability would remain negligi- 
bly small. 

5. Electromotive force for AMRI 

Stratorotational instability under the influence of a toroidal field 
is tempting to apply the concept of the mean-field electrodynam- 
ics in turbulent fields. The nonaxisymmetric components of flow 
and field can be used as fluctuations while the axisymmetric 
components are considered as the mean quantities. Simply the 
averaging procedure is the integration over the azimuth <p. 

It is standard to express the turbulence-induced electromo- 
tive force as 



£ = (u' X B') = a o (B) - 77 T curl<5> 



(21) 



with the alpha-tensor a and the (scalar) eddy diffusivity tj-y. 
In cylindric geometry the mean magnetic field (B) has only 
a ^-component and the mean current curl(B) only has a z- 
component. Hence, £^ = a^{B^) and & z = —7]tC\xt\{B) = 0. 
The latter relation only holds for AMRI where the mean mag- 
netic field is current-free. 

In the present paper we only consider the alpha-effect in the 
frame of a linear theory where all functions are free to one and 
the same (complex) arbitrary parameter. This is only possible 
if the second-order quantities such as £^ are normalized with a 
second-order quantity. In the considered case it makes sense to 
form the correlation function 



r 



<h' x B') 



MAX( y/(u' 2 )(B' 2 }) 



(22) 



where all the correlations along the radius R are normalized with 
one and the same parameter. Hence, this function does no longer 
contain the arbitrary factor of the eigenfunctions and by defini- 
tion it is smaller than unity. The main questions are the sign of 
this quantity and the influence of the density stratification. Note 
that the has been given as positive in the model setup where 
rotation axis and density gradient are antiparallel. The latter is 
quite characteristic for the situation at the poles of rotating stars 
or disks so that the old question whether an alpha-effect at the 
poles does exist or not is here concerned. 

One can easily show that second-order correlations of quan- 
tities running with exp(i(fe + m<f>)) after integration over (f> do 
not depend on the coordinate z. The term kz only fixes the phase 
where the integration starts. 

The correlation function f" can be estimated for fast-rotating 
magnetoconvection in the high conductivity limit a =s ~<J(u' 2 } 
and with (B' 2 } * {qj I rf){B) 2 after Krause & Riidler (1980) so 
that 



f < 



(23) 



For stellar convection this is a small number. Hence, the numer- 
ical constraints for the exact calculations of the eigenfunctions 
are very high. One must keep this finding in mind when dis- 
cussing the numerical results given in Fig. [12] for an example 
without density stratification and another one with density strat- 
ification. For simplicity the magnetic Prandtl number is put to 
unity. Note the influence of the density stratification. The ex- 
amples are taken from the Fig. [4] for Ha = 100. The resulting 



G. Riidiger & D.A. Shalybkov: Stratorotational instability in MHD Taylor-Couette flow 



9 




-0.2 



Ha=100 
-Ha=200 




Ha=400 



RADIUS 



RADIUS 




RADIUS 

Fig. 12. The correlation function ( l22l > for AMRI with conduct- 
ing cylinders and with fj = 0.5, fin = 0.35 (pseudo-Kepler), 
fiB = 0.5, Pm = 1. Density stratification is zero (top) and finite 
(bottom, Fr = 0.5). The Reynolds numbers Re = 141 (top) and 
Re = 341 (bottom) were taken from Fig. [4] 



correlations for N - are much more antisymmetric than the 
values for Fr = 0.5. Averaged over the radius the f a vanishes for 
N = but not for the presence of a density stratification. 

This cannot be a boundary effect. A boundary effect is im- 
plausible as the preferred radial direction is perpendicular to the 
rotation axis for cylinders which mimics the equator in rotating 
spheres and there the c-effect must vanish. Indeed, the antisym- 
metry is reduced if a density stratification is allowed (Fig. [12J 
bottom). Then the density gradient is now the preferred direc- 
tion in the system parallel to the rotation allowing the formation 
of large-scale helicity (as in rotating spheres at the poles). We 
are thus tempted to assume that the presented calculations for 
the presence of a density gradient indeed have demonstrated for 
the existence of an a-effect for the magnetic SRI. 

Figure [13]presents similar results but for the small magnetic 
Prandtl number Pm = 0.01. The correlation functions do pro- 
vide a dominance of the positive contributions but for higher and 
higher magnetic fields the difference becomes smaller. It seems 
that the numerical results do here reflect a magnetic quenching 
of the a-=-effect. 



Fig. 13. The correlation function d22l) for fj = 0.5, fin = 0.35 
(pseudo-Kepler), fig = 0.5, Pm = 0.01. Fr = 0.5). The Reynolds 
numbers are Re = 531 for Ha = 100 and Re = 999 for Ha 
200 and Re = 2137 for Ha = 400. Note that the alpha-effect is 
quenched for stronger magnetic fields. 

6. Conclusions 

The stability of the dissipative Taylor-Couette flow under the 
joint influence of a stable vertical density stratification and an 
azimuthal magnetic field is considered. The problem is of inter- 
est for future laboratory experiments but also within the frame 
of accretion disk physics. The Kepler rotation generates strong 
toroidal magnetic fields dominating the poloidal components. 
The standard MRJ which works with only axial fields may be 
of minor relevance for the stability of the Kepler rotation law 
compared with the azimuthal MRJ in connection with the den- 
sity stratification and Tayler instability. Mainly nonaxisymmet- 
ric 'kink' modes (m — 1) are considered but there are also exam- 
ples where the axisymmetric modes are the most unstable ones. 

We started with a discussion of the SRI without magnetic 
fields. For a flat rotation law a stratification value Fr exists for 
which the critical Reynolds number of the rotation has a min- 
imum (see Fig. The SRI is basically stabilized by both too 
weak or too strong stratification. The instability only appears if 
the characteristic buoyancy time approaches the rotation period. 
Also if the rotation is too flat the instability disappears. 

The limiting ratio fin strongly depends on the gap width (see 
Fig. [3). For small gaps rotation laws with fin > fj even prove 
to be unstable while for wide gaps the condition fin < H results 
for exciting SRI. Our previous finding that fin ~ f) limits the SRI 
is reproduced for medium-sized gaps. In all cases, however, the 
quasi-Kepler rotation proves to be unstable. New experiments 
with different gap sizes and larger Reynolds numbers could ver- 
ify these results. 

Figures [4] and [5] yield the basic results for SRI subject to 
toroidal fields. The magnetic field is assumed as current-free in 
the fluid between the cylinders, i.e. B$ oc l/R excluding Tayler 
instability. Without density stratification no Rayleigh instabil- 
ity exists for quasi-Kepler rotation but nonaxisymmetric modes 
with m = 1 are unstable as a result of the interaction of differen- 
tial rotation and magnetic field ('AMRI'). 

The magnetic influence strongly depends on the magnetic 
Prandtl number. The instability needs higher Reynolds number 
for Pm < 1 and it needs lower Reynolds number for Pm > 1 . 
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After our experiences with MHD instabilities this is not a sur- 
prise. It is a surprise, however, that always the magnetic influ- 
ence is only weak. Up to Hartmann numbers of Ha 100 only 
a magnetic-induced factor of two plays a rolfl Hence, our con- 
clusion is that the SRI survives for rather high magnetic fields. 
For large Pm it is even supported by the toroidal magnetic field. 

The combination of rotation, density stratification and mag- 
netic field leads to complex results. However, a basic observation 
is that the critical Reynolds numbers, if existing, above which the 
flow becomes unstable without a magnetic field, are increased by 
the stable density stratification. In contrast, the critical Hartmann 
numbers, if existing by the Tayler instability, above which the 
field becomes unstable without a rotation, do not depend on the 
stratification. Different transition are possible between the two 
limits. 

Generally speaking, the density stratification 'balloons' the 
stability region and in this sense it stabilizes the flow. For slow 
rotation the maximal stable magnetic field exceeds the critical 
magnetic field without rotation while for faster rotation the max- 
imal stable magnetic field is smaller than this critical value. Even 
rather slow values of the Reynolds numbers lead to a stabiliza- 
tion of those fields which are unstable for Re = 0. The effect is 
strong for Pm = 1 but it becomes smaller for decreasing mag- 
netic Prandtl numbers. 

For steep radial profiles of the magnetic field (i.e. strong ax- 
ial currents) and magnetic Prandtl numbers Pm > 1 one also finds 
the magnetic field destabilizing the Rayleigh instability, i.e. the 
critical Reynolds numbers with magnetic field are lower than the 
Reynolds numbers without magnetic field. This magnetic desta- 
bilization only exists for not too small magnetic Prandtl num- 
bers. It exists for both uniform and density-stratified fluids (Fig. 
ED- 

Finally, the ^-component of the electromotive force repre- 
senting the a-effect of the mean-field dynamo theory has been 
computed. The computations require an extreme degree of ac- 
curacy. The results demonstrate the importance of the density 
stratification for the existence of the a-effect. Without density 
stratification the correlations are vanishing in the radial average. 
With an axial density stratification the calculations model the 
polar region of a rotating sphere or disk. This interpretation ac- 
cepted we found the at the northern pole as positive (Fig.fTZli. 
Again, the basic ingredient of this a-effect model is the density 
stratification. 
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2 Ha =100 for gallium corresponds to B as 2200/Wo [Gauss] with^o 
in cm 



